Deterministic scale-free networks created in a recursive manner 
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In a recursive way and by including a parameter, we introduce a family of deterministic scale-free 
networks. The resulting networks exhibit small-world effects. We calculate the exact results for 
the degree exponent, the clustering coefficient and the diameter. The major points of our results 
indicate: the degree exponent can be adjusted; the clustering coefficient of each individual vertex is 
inversely proportional to its degree and the average clustering coefficient of all vertices approaches 
to a nonzero value in the infinite network order; and the diameter grows logarithmically with the 
number of network vertices. 

PACS numbers: 02.10.Ox, 89.75.Hc, 89.75.Da, 89.20.Hh 



I. INTRODUCTION 

Since the two seminal papers by Watts and Strogatz 
on small-world networks [lj and Barabasi and Albert on 
scale- free networks Q , we have witnessed a considerable 
efforts devoted to characterization and understanding of 
complex networks, which describe many systems in na- 
ture and society. In the past few years, there has been a 
substantial amount of interest in network structure and 
function from a wide circle of researchers [E 0, IE IE • 
One particular question that has attracted an exceptional 
amount of attention concerns the structure of networks 
that are evolving over time. A number of network models 
have been proposed, which convincingly reproduce some 
or all features of real- life systems [E 0,1E IE Q • 

Most of the precious models are stochastic [E El IE IE 0- 
However, because of their advantages, deterministic net- 
works have also received much attentio n pi IE fToL ITU 

EEEllIEIiEHEEHH EElH. First, 

the method of generating deterministic networks makes it 
easier to gain a visual understanding of how networks are 
shaped, and how do different vertices relate to each other; 
moreover, deterministic networks allow to compute ana- 
lytically their properties: degree distribution, clustering 
coefficient, average path length, diameter, betweenness, 
modularity and adjacency matrix whose eigenvalue spec- 
trum characterizes the topology H H HE [lj EH III 

he he m in m he in he m 

The first model for deterministic scare-free networks 
was proposed by Barabasi et al. in Ref. Q and was 
intensively studied in Ref. Q. Another el egan t model, 
called pseudofractal scale-free web (PSW) [l0|. was in- 
troduced by Dorogovtsev and Mendes, and was extended 
by Cornelias et al. in Ref. ^3 • Based on a similar idea of 
PSW, Jung et al. presented a class of recursive trees [l2j . 
Additionally, in order to discuss modular ity. Ravasz et al. 
proposed a hierarchical network model (l3l IT^j , the ex- 
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act scaling properties and extensive study of which were 
reported in Refs. HE an d EE> respectively. Recently, In 
relation to the problem of Apollonian space-filing pack- 
ing, Andrade et al. introduced Apollonian networks Qj] 
which were also proposed by Doye and Massen in Ref. Hal 
and have been intensively investigated 0, I^E I2H |22| . 
Except for above models, deterministic networks can be 
created by various techniques: modification of some reg- 
ular graphs l23l. addition and product of graphs |24j| . 
edge iterations J25j and other mathematical methods as 
in Ref. Hi 

As mentioned by Barabasi et al., it would be of ma- 
jor theoretical interest to construct deterministic models 
that lead to scale-free networks |8(. In this paper, we do 
an extensive study on pseudofractal scale-free web ^3 
and the recursive graphs (RG) ^3- I n a simple recur- 
sive way we propose a general model for PSW and RG 
by including a parameter, with PSW and RG as partic- 
ular cases of the present model. The deterministic con- 
struction of our model enables one to obtain the analytic 
solutions for some main structure properties: degree dis- 
tribution, clustering coefficient and diameter. By adjust- 
ing the parameter, we can obtain a variety of scale-free 
networks. 



II. THE NETWORK CONSTRUCTION 

Before introducing our networks we give the following 
definitions on a graph. The term size refers to the num- 
ber of edges in a graph. The number of vertices in a 
graph is called its order. When two vertices of a graph 
are connected by an edge, these vertices are said to be 
adjacent, and the edge is said to join them. A complete 
graph is a graph in which all vertices are adjacent to one 
another. Thus, in a complete graph, every possible edge 
is present. The complete graph with q vertices is denoted 
as K q (also referred in the literature as q- clique; see pTi). 
Two graphs are isomorphic when the vertices of one can 
be relabeled to match the vertices of the other in a way 
that preserves adjacency. So all g-cliques are isomorphic 
to one another. 
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FIG. 1: Scheme of the growth of the deterministic network 
for the case of q — 2 and m = 2. Only the first three steps 
are shown. 
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respectively. So for large t, The average degree k t = - J jy i£ 
is approximately 2q. 



III. CHARACTERISTICS OF THE NETWORKS 
A. Degree distribution 

When a new vertex i is added to the graph at step ti, 
it has degree q and forms q new g-cliques. Let n q (i, t) be 
the total number of g-cliques at step t that will created 
new vertices connected to the vertex i at step t + 1. So 
at step ti, n q (i, ti) = q. By construction, we can see that 
in the subsequent steps each new neighbor of i generated 
q — 1 new g-cliques with i as one vertex of them. Let 
ki(t) be the degree of i at step t. It is not difficult to find 
following relations for t > ti + 1: 



Afcj(t) = k^t) - fc f (t - 1) = mn q (i, t - 1) (3) 



and 



The network is constructed in a recursive way. We 
denote the network after t steps by R(q,t), q > 2,t > 
(see Fig. Then the network at step i is constructed 
as follows: For t = 0, i?(g, 0) is a complete graph A" 9+ i 
(or (g + l)-clique) consist of q+ 1 g-cliques), and i?(g, 0) 
has g + 1 vertices and g(g + l)/2 edges. For t > 1, i?(g, t) 
is obtained from i?(g, t — 1) by adding m new vertices for 
each of its existing subgraphs isomorphic to a g-clique, 
and each new vertex is connected to all the vertices of this 
subgraph. In the special case m = 1, it is reduced to the 
network described in Ref. which is a generalization 
of pseudofractal scale-free web [Tof . 

Let n v (t) and n e (t) be the number of vertices and 
edges created at step t, respectively. Denote K q _ t as 
the total number of g-cliques in the whole network at 
step t. Note that the addition of each new vertex leads 
to q new g-cliques and g new edges. By construc- 



tion, we have n e (t) — qn v (t), n v {t) 



mK, 



K q ,t = 
K q ,t = 



q,t-l 



and 



Kq,t-i + qn v (t). Thus one can easily obtain 
(mq + l)K q , t -i = (q + l)(mq + if (t > 0), 
n v (t) = m(q + l)(mq + l)*" 1 (t > 0) and n e (t) = 
mq(q + l)(mq + (t > 0). So the number of net- 

work vertices increases exponentially with time. From 
above results, we can easily compute the size and order 
of the networks. The total number of vertices N t and 
edges \E\ t present at step t is 

AT us (g + l)[(mg + 1^+9-1] 

N t = 2_^n v (ti) = (1) 



ti=0 



and 



\E\ t = ^ n e (U) = (q + l)(mq + 1)' + 2) ,(2) 



n q (i,t) = n q (i,t-l) + (q- l)Aki{t). 



(4) 



From the above two equations, we can derive n q (i, ti) — 
[m(q — 1) + l]n q (i, ti — 1). Considering n q (i,ti) = q, 
we obtain n q (i,t) — q[m(q — 1) + l] t_4i and Afcj(t) = 
mq[m(q — 1) + l] t_ * i_1 . Then the degree ki(t) of vertex 
i at time t is 



h(t) = ki(U)+ Ak *(th) 
t h =n+i 

q[m{q - 1) + I]*'*' + q 2 - 2g 
9-1 



(5) 



Since the degree of each vertex has been obtained ex- 
plicitly as in Eq. ifHjl. we can get the degree distribu- 
tion via its cumulative distribution, i.e. P C um(k) = 
J2k'>k N ( k '> t )/ N t ~ fcl ~ 7 ; where N(k',t) denotes the 
number of vertices with degree kl . The analytic compu- 
tation details are given as follows. For a degree k 



q[m(q- 1) + l]'" 3 + q 2 - 2q 



(6) 



there are n v (j) = m(q + l)(gm + vertices with this 
exact degree, all of which were born at step j. All vertices 
with birth time at j or earlier have this and a higher 
degree. So we have 



J2N(k',t)=J2n v (s) 



(g + l)[(mg + lp+g-l] 



k'>k 



s=0 



q 



As the total number of vertices at step t is given in Eq. (JIJ 
we have 



t;=0 



g[m(g-l) + l]«-J+g 2 -2g 
9-1 



1-7 (q+l)l(rnq+iy+q-l] 
_ 9 

(q+l)l( m q+iy+q-l] 
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Therefore, for large t we obtain 



[[m(?-l) + l]*-*] X 7 = H + l) i_i 



and 



7 f» 1 + 



hi(mq + 1) 
ln[m(q - 1) + 1] ' 



(7) 



(8) 



For the particular case of to = 1, Eq. (JSJ) recovers the 
results previously obtained in Ref. [ill ]. 



B. Clustering coefficient 

The clustering coefficient [l| Ci of vertex i is defined 
as the ratio between the number of edges that actually 
exist among the ki neighbors of vertex i and its maximum 
possible value, ki(ki — l)/2, i.e., Ci = 2ei/ki(ki — 1). The 
clustering coefficient of the whole network is the average 
of C[s over all vertices in the graph. 

For our networks, the analytical expressions for cluster- 
ing coefficient C{k) of the individual vertex with degree 
k can be derived exactly. When a vertex is created it is 
connected to all the vertices of a q-clique whose vertices 
are completely interconnected. Its degree and cluster- 
ing coefficient are q and 1, respectively. In the following 
steps, if its degree increases one by a newly created ver- 
tex connecting to it, then there must be q — 1 existing 
neighbors of it attaching to the new vertex at the same 
time. Thus for a vertex of degree k, we have 



<?(<?- 



C{k) 



il + ( g -l)(*-g) 2 (<Z-l)(fc-§) 



fc(fc-l) 



k(k - 1) 



(9) 



which depends on degree k and q. For k 3> q, the C(k) is 
inversely proportional to degree. The scaling C(k) ~ k~ l 
has been found for some network models [ToL ITU IbH . ITR 
IrH ITsl Pill I20I I21I , and has also observed in several real- 
life networks [l4| . 

Using Eq. , we can obtain the clustering Ct of the 
networks at step t: 



C t = 



t 

E 

r=0 



2(q-l)(D r 



f)n„(r) 



L» r (D r - 1) 



(10) 



where the sum is the total of clustering coefficient for all 
vertices and D r = Umiq-^+i^+q 2 -^ shown by Eq _ p 

is the degree of the vertices created at step r. 

It can be easily proved that for arbitrary fixed to, Ct in- 
creases with q, and that for arbitrary fixed q, Ct increases 
with to. In the infinite network order limit (N t — > 00), 
Eq. 110(1 converges to a nonzero value C. When q = 2, 
for to = 1, 2, 3 and 4, C equal to 0.8000, 0.8571 0.8889 
and 0.9091, respectively. When to = 2, for q = 2, 3, 4 and 
5, C are 0.8571, 0.9100, 0.9348 and 0.9490, respectively. 
Therefore, the clustering coefficient of our networks is 
very high. Moreover, similarly to the degree exponent 7, 
clustering coefficient C is determined by q and to. Fig. 
121 shows the dependence of C on q and to. 




FIG. 2: The dependence relation of C on q and m. 



C. Diameter 

The diameter of a network characterizes the maximum 
communication delay in the network and is defined as the 
maximum of shortest path between all pairs of vertices. 
In what follows, the notations \x] and [^J express the 
integers obtained by rounding x to the nearest integers 
towards infinity and minus infinity, respectively. Now we 
compute the diameter of R(q,t), denoted Diam(R(q,t)) 
for q > 3 (q = 2 is a particular case that is treated 
separately) : 

Step 0. The diameter is 1. 

Steps 1 to [^Tp-] — 1- In this case, the diameter is 
2, since any new vertex is by construction connected 
to a g-clique forming a (q + l)-clique, and since any 
(q + l)-clique during those steps contains at least f 2 ^-] 
(q even) or [^^-]+l {q odd) vertices from the initial 
(q + l)-clique R(q, 0) obtained after step 0. Hence, any 
two newly added vertices u and v will be connected 
respectively to sets S u and S v , with S u C V(R(q 7 0)) 
and S v C V(R(q,0)), where V(R(q, 0)) is the vertex set 
of R(q,0); however, since \S U \ > [^rl (<Z even) and 
\S V \ > (q odd), where \S\ denotes the num- 

ber of elements in set S, we conclude that S u n S v ^ 0, 
and thus the diameter is 2. 

Steps to q. In any of those steps, some newly 

added vertices might not share a neighbor in the original 
(q + l)-clique i?(g,0); however, any newly added vertex 
is connected to at least one vertex of the initial (q + 1)- 
clique R(q,0). Thus, the diameter is equal to 3. 

Further steps. In order to simplify the analysis, we 
first note that it is unnecessary to look at all the ver- 
tices in the graph in order to find the diameter. In other 
words, some vertices added at a given step can be ig- 
nored, because they do not increase the diameter from 
the previous step. These vertices are those that connect 
to vertices that already formed a (q + l)-clique in one or 
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several of the previous steps: indeed, for these vertices we 
know that a similar construction has been done in pre- 
vious steps, so we can ignore them for the computation 
of the diameter. Let us call "outer" vertices the vertices 
which arc connected to a g-clique that did not exist pre- 
viously. Clearly, at each step, the diameter depends on 
the distances between outer vertices. 

Now, at any step t > q+l, we note that an outer vertex 
cannot be connected with two or more vertices that were 
created during the same step < t' < t — 1. Moreover, 
by construction no two vertices that were created during 
a given step are neighbors, thus they cannot be part of 
the same g-clique. Therefore, for any step t > q+l, some 
outer vertices are connected with vertices that appeared 
at pairwise different steps. Thus, if Vt denotes an outer 
vertex that was created at step t, then vt is connected to 
vertices UjS, 1 < i < t — 1, where all the is are pairwise 
distinct. We conclude that vt is necessarily connected to 
a vertex that was created at a step to < t — q. If we 
repeat this argument, then we obtain an upper bound on 
the distance from vt to the initial (q + l)-clique R(q,0). 
Let t = aq + p, where 1 < p < q. Then, we see that 
Vt is at distance at most a + 1 from a vertex in R(q,t). 
Hence any two vertices Vt and u>t in R(q, 0) lie at distance 
at most 2{a + 1) + 1 ; however, depending on p, this 
distance can be reduced by 1, since when p < [2i-] — 1, 
we know that two vertices created at step p share at least 
a neighbor in R(q,0). Thus, when 1 < p < [2+1] - 1, 
Diam(R(q,t)) < 2(a + 1), while when [2±1] < p < q, 
Diam(R(q,t)) < 2(a + 1) + 1. One can see that these 
bounds can be reached by pairs of outer vertices created 
at step t. More precisely, those two vertices vt and Wt 
share the property that they are connected to q vertices 
that appeared respectively at steps t — l,t — 2, , . A — q. 

Based on the above arguments, one can easily see that 
for t > q, the diameter increases by 2 every q steps. More 
precisely, we have the following result, for any q > 3 and 
t > 1 (when t = 0, the diameter is clearly equal to 1): 

Diam(K(q,t)) = 2([—\ + 1) + f(q, t), 

q 

where f(q,t) = if t - [^\q < [2±1] - 1, and 1 oth- 
erwise. When t gets large, Diam(R(q,t)) ~ y, while 
N t ~ (mq + 1)*, thus the diameter grows logarithmically 



with the number of vertices. 

For q = 2, the argument is similar: when t = (resp. 
t = 1), the diameter is equal to 1 (resp. 2). Now consider 
two outer vertices created at step t > 2, say vt and Wt- 
Then vt is connected to two vertices, and one of them 
must have been created before or during step t — 2. We 
repeat this argument, and we end up with 2 cases: (1) 
t = 2m is even. Then, if we make m "jumps" , from Vt 
we arrive in i?(2,0), in which we can reach any wt by 
using an edge of i?(2, 0) and making m jumps to Wt in 
a similar way. Thus Diam(R(2,2m)) < 2m + 1. (2) 
t = 2m + 1 is odd. In this case we can stop after m 
jumps at i?(2,l), for which we now that the diameter 
is 2, and make m jumps in a similar way to reach u> t . 
Thus Diam(R(2,2m + 1)) < 2m + 2. As previously, it 
easily seen that the bound can be reached for some pairs 
of vertices. 

Hence, formally, Diam(R(2,t)) = t + 1 for any t > 0. 
In this particular case, the network order Nt ~ (2m+ 1)' , 
thus Diam(R(2,t)) also increases logarithmically with 
network order. 



IV. CONCLUSION AND DISCUSSION 

In summary, we have proposed and studied a network 
model constructed in a recursive fashion. At each time 
step, each already existing q-clique generates m new ver- 
tices. The process results in a serial of networks with two 
published papers [Tol . as special cases of them. We 
have obtained the analytical results for degree exponent, 
clustering coefficient and diameter. The degree exponent 
and the clustering coefficient may be adjusted to various 
values by tuning the parameter m. Therefore, they may 
perform well in mimicking a variety of scale-free networks 
in real world. Moreover, in a similar way, one can easily 
consider other variations, for example, at each step not 
all cliques of a network but only some of them generate 
vertices, which also allow a rich structure and flexibility 
in the control of degree exponent, clustering coefficient 
and other properties. Details of the analytical solution 
for these variants will be addressed elsewhere. 

This research was supported by the National Natural 
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